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Abstract 

We consider real random symmetric N x N matrices H of the band- 
type form with characteristic length b. The matrix entries H(x,y),x < y 
are independent Gaussian random variables and have the variance propor- 
tional to where u(t) vanishes at infinity. We study the resolvent 
G(z) = (H — z)~ , Imz 7^ in the limit 1< i C JV and obtain explicit 
expression 5(2:1,2:2) for the leading term of the first correlation function 
of the normalized trace (G(z)> = N' 1 Tr G(z). 

We examine S(Xi +i0, A2 — i0) on the local scale Ai — A2 = and show 
that its asymptotic behavior is determined by the rate of decay of u(t). 
In particular, if u{t) decays exponentially, then S(r) ~ —Cb 2 N~ 1 r~ 3 ^ 2 . 
This expression is universal in the sense that the particular form of u de- 
termines the value of C > only. Our results agree with those detected in 
both numerical and theoretical physics studies of spectra of band random 
matrices. 

1 Problem, motivation and results 

Random matrices play an important role in various fields of mathematics and 
physics. The eigenvalue distribution of large matrices was initially considered 
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by E.Wigncr to model the statistical properties of energy spectrum of heavy 
nuclei (see e.g. the collection of early papers Q). Further investigations have 
led to numerous applications of random matrices of infinite dimensions in such 
branches of theoretical physics as statistical mechanics of disordered spin sys- 
tems, solid state physics, quantum chaos theory, quantum field theory and others 
(see monographs and reviews ||, |l^, [l6|, Q). In mathematics, the spectral the- 
ory of random matrices has revealed deep links with the orthogonal polynomi- 
als, integrable systems, representation theory, combinatorics, non-commutative 
probability theory and other theories |j, 11, [32, 35 1. 

In present paper we deal with the family of real symmetric random matrices 
that can be referred to as the band-type one. In the simplest case the matrices 
have zeros outside of a band around the principal diagonal. Inside of this band 
they are assumed to be jointly independent random variables. The limiting 
transition considered is that the band width b increases at the same time as the 
dimension of the matrix n does. 

There is a large number of papers devoted to the use of random matrices of 
this type in models of quantum chaotical systems (see, e.g. Q and references 
therein). In these studies, one of the central topic is related with the transition 
between fully developed chaos and complete integrability. The crucial observa- 
tion made numerically |j) and then supported in the welth of theoretical physics 
papers (see, for example |33|]) is that the ratio b 2 /n is the critical one for 
the corresponding transition in spectral properties of band random matrices. 

On the rigorous level, the eigenvalue distribution of if( ra > b ) has been studied 
[[|, ||, |2(|. It is shown that the limiting eigenvalue distribution exists, is non- 
random and depends on the parameter a = linin^oo b/n. However, the role 
of the ratio b 2 /n has not been revealed there. Recently, a series of papers 
appeared where the band random matrices are studied in the context of the 
non-commutative probability theory jl7| These studies also deal with the 
limit n, b — > oo such that a > 0. 

In present paper we are concentrated on the case of a = represented by 
the limit 

1 < b < n 

and study the first correlation function of the resolvent of band random matrices. 
We show that the ratio /3 = linin^oo b 2 /n naturally arises when one considers 
the leading term of this correlation function on the local scale. This can be 
regarded as the support of the conjecture that the local properties of spectra of 
band random matrices depend on the value of 0. 

Let us describe our results in more details. We consider the ensemble 
{H( n ' b '} of random N x N matrices, N = 2n + 1 whose entries H^ n ' h \x,y) 
have the variance proportional to m(^t^)j where u(t) > vanishes at infinity. 

We consider the resolvent G^ n ' b ^ (z) = (iJ (™> b ) — z) 1 and study the asymptotic 
expansion of the correlation function 

C nt b(zi,Z2) = ~E 1 fn,b{zi)fn,b{z2) ~ E/ n)6 (2i)E/ n)fe (^ 2 ), 
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where we denoted f n ,b(z) = N^ 1 Tr G^ n,b '(z). Keeping z% far from the real 
axis, we consider the leading term S(zi,Z2) of this expansion and find explicit 
expression for it. This term S(ri + i0,r2 — iO) regarded on the local scale 
n — T2 = r/N exhibits different behavior depending on the rate of decay of the 
profile function u(t). 

Our main conclusion is that if u(t) ~ \t\ v as t — > oo, then the value v = 2 
separates two major cases. If v G (1, 2), then the limit of 5(r) depends on v. If 
f G (2, +oo), then 

]V6 5W — COnSt -^'FF (1 + 0(1)) - 

These results are in agreement with those predicted in theoretical physics stud- 
ies. In particular, the last expression for S coincides with the Altshulcr-Shklovski 
asymptotics of the spectral correlation function (see e.g. p7[). 

The paper is organized as follows. In Section 2 we determine the family of 
ensembles and present several already known results that will be needed. In 
Section 3 we formulate our main propositions and describe the scheme of their 
proofs. To illustrate this scheme, we present a short proof of the Wigner semi- 
circle law for Gaussian Orthogonal Ensemble of random matrices. In Section 4 
we study the correlation function C n ,b(zi, z%) and obtain the explicit expression 
S(z\, zz) for its leading term. In Section 5 we study the self-averaging property 
of G(z) and prove auxiliary facts used in Section 4. Expressions derived in Sec- 
tion 4 are analyzed in Section 6, where the asymptotic behavior of 22) is 
studied. In Section 7 we give a summary of our observations. 



2 Band Random Matrices and Wigner Law 
2.1 The ensemble 

Let us consider the family A — {a(x, y), x < y, x, y € Z} of jointly independent 
random variables determined on the same probability space. We assume that 
they have joint Gaussian (normal) distribution with properties 

Ea(x,y) = 0, E [a (a;, y)] 2 = v(l + S xy ), (2.1) 
where we denote by S the Kronecker symbol; 



S X y — 



0, if x ^ y, 

1, if x = y. 



Here and below E denotes the mathematical expectation with respect to the 
measure generated by the family A. 

Let u(t),t S R be a piece-wise continuous function u(t) = u(—t) > satis- 
fying conditions 

sup \u(t)\ — u < 00 (2.2) 
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and 

u(t)dt = 1. (2.3) 
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For simplicity, we assume u(t) to be monotone for t > 0. 

Given real parameter b > 0, we introduce an infinite matrix U^' 

£/W( x , y ) = i u fcj^, x,yeZ. 

and determine the ensemble as the family of real symmetric matrices 

of the form 

H^(x,y) = a(x,y)yJuW(x,y), x<y, \x\, \y\ < n, (2.4) 

where b < N, N = 2n + 1 and the square root is assumed to be positive. 

Let us note that the matrix (2.1) has the really band form when is 
constructed with the help of a function u having a finite support, say 



u(t) 



l, ifte(-±,±), 

0, otherwise. 



In this case the band width is less than or equal to 26 + 1. If b = N, then 
matrices (2.4) coincide with those belonging to Gaussian Orthogonal Ensemble 
(GOE) (25|. This random matrix ensemble is determined as the family {An} 
of real symmetric matrices 

A N (x,y) = —=a(x,y), x,y=l,...,N, (2.5) 
v W 



with {a(x,y)} belonging to A (2.1). GOE together with its Hermitian and 
quaternion versions plays the fundamental role in the spectral theory of random 
matrices. 

Random symmetric matrices (2.5) with independent arbitrary distributed 
random variables a(x, y) satisfying (2.1) is referred to as the Wigner ensemble of 
random matrices. This random matrix ensemble considered first by E. Wigner 
[ j36| is extensively studied in a series of papers (see e.g. H|| and references 
therein). In particular, in paper ]22] the resolvent technique is developed to 
study the spectral properties of the Wigner ensemble. Actually, we follow a 
version of this technique, but restrict ourself with more simple case of Gaussian 
random variables. More general case of arbitrary distributed random variables 
would make the computations much more cumbersome. Let us repeat that the 
main task of this paper is to study the role of the ratio between b and N with 
respect to the spectral properties of random matrices. 

Finally, it should be noted that we restrict ourself with the ensemble of real 
symmetric matrices for the sake of simplicity also. All results can be obtained 
by using essentially the same technique for the Hermitian analogue of H^ n ' b \ 
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2.2 Limiting eigenvalue distribution 

Eigenvalue distribution of matrices if ( n ' & ) is described by the normalized eigen- 
value counting function 

ct(A; /f ( "< b) ) = #{\ { "' b) < \}N-\ (2.6) 

where A^' b) < . . . < A^ 6) are eigenvalues of H^ n ' b ' . We denote by fn,b(z), z <G 
C the Stieltjes transform of the measure given by (2.6); 

/„,&(*) = r Imz^O. (2.7) 

Given a Stieltjes transform f(z), one can restore corresponding measure da(X) 
with the help of the inversion formula (see e.g. 

The limiting behavior of (2.7) as n, b — ► oo was studied in a series of papers 



23, E6|. It was proved in [E6l that f n ,b(z) converges as ^ oo in 



probability to a nonrandom function that depends on the ratio a = ]hnb/N; 

p- lim f n ,b(z) = w a (z). (2.8) 

n,b — >oo 

In particular, if a = 0, then the function wq(z) = w(z) satisfies equation 

w(z) = X — - (2.9) 

— z — vw(z) 

The solution of this equation is unique in the class of functions satisfying con- 
dition 

Imw(z)Imz > 

and can be represented in the form w(z) — J (A — z)~ 1 da w (X), where cr^A) is 
the famous semicircle (or Wigner) distribution [^6| with the density 

r\\ i { \\ 1 f V4w - A 2 , if |A| 2 < 4u, /oin\ 
^(A)=^(A) = — | Q) if j A j 2 > 4w ; ( 2 -!0) 

This density was obtained first by E. Wigner |36| ] for eigenvalues of random 
matrices of the "full" form (2.5) and can be also obtained as the limit (2.8) with 
a = 1 

^(A) = lim a{\;A N ) (2.11) 

N— *oo 

Thus, one gets the same eigenvalue distribution in the opposite limiting 
transitions of narrow a — and wide a = 1 band widths. It is known that 
in the intermediate regime < a < 1 the limiting distribution differs from 
the semicircle (2.11) 26 . In present paper we concentrate ourself on the most 
interesting case a — 0. 
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In present paper we always consider the case of 1 <C b <C n. As we have 
noted, the paper is aimed to detect the role of the parameter [3 = liniAr^oo b 2 /N. 
To avoid technical problems, we restrict ourself with the range 

b = n x 1/3<x<1- (2.12) 

We are convinced that our results are valid on the whole range < \ < 1. 

3 Main Propositions and Scheme of the Proof 

The resolvent 

G^ b) (z) = - ziy 1 , Imz ^ 

is widely exploited in the spectral theory of operators. Its normalized trace 
(G^ n ' b \z)) coincides with the Stieltjes transform f n ,b{z) (2.7); 

The results of this section are related with the asymptotic behavior of 
(G("> b ) (z)) in the limit (2.12), with z G A,,, 

A, = {z e C : |Imz| > 77} with 77 = + 1. (3.1) 
3.1 Main technical results 

Our first statement concerns the pointwise convergence of the diagonal entries 
Qi n fi){x,x; z), \x\ < n of the resolvent. Let us determine the set 

B L = B L (n,b) = {x e Z : \x\ < n- bL} . (3.2) 

Theorem 3.1 

Given e > 0, there exists a natural L such that 

sup \G {n > b) (x,x;z)-w(z)\ <e, Vz e A„ (3.3) 

xeB L 

for sufficiently large b,n. 

The result of Theorem 3.1 is interesting by itself. We shall use it hardly in 
the proof of the following statement concerning the correlation function 

C n ,b(z lt Z2) = E (G(*i)) (G(z 2 )) - E (G(*i)) E (G(z 2 )) . (3.4) 
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Theorem 3.2. 

If Zi G A^, i = 1,2, then in the limit n, b — > oo (2.12) 



C nth {z u z 2 ) = ■^rS(z 1 ,z 2 ) + o(^). (3.5) 



The explicit term of S{z\ 1 z 2 ) is given by relation 

2v 

(1 — vwf) (1 — vw\) 
where Wj = w(zj), j = 1, 2 and 2:2) is given by the formula 

i>k;iu;2«f(p)] 

where we denote by uf(p) the Fourier transform of u 



S(zi, z 2 ) = 2 7- T -Q(zi, z 2 ), (3.6) 



Q{Zl,Z2) = — - _ 2 dp, 



{^(p) = / u(i)e ipt di. 



R 



It should be noted that in the case of GOE (2.6) relation (3.5) is valid with 
b replaced by N and expression (3.6) takes the following form (see e.g. |l4|, p0[ ) 

Agoe(Z1 ' Z2) (l-vwl)(l-v W i 2 ) {1 - VWlW2 f (3J) 

Let us briefly explain why (3.6) and (3.7) lead to different asymptotic expressions 
on the local scale determined as 

^A+^+iO, zf^A + ^-iO (3.8) 

with A 6 supp da w (2.10). It follows from equality (2.9) that 

2 

'JlW 2 I wi — w 2 ' 



[1 — 0W1W2] 2 \ z\ — z 2 



(3.9) 



This expression tends to infinity in the limit (3.8) and vw{z\)w{z 2 ) — > 1 as well. 
But after dividing by iV 2 , one obtains from (3.7) and (3.9) that 

-^S GOE ( Zl ,z 2 ) = — L^(l + (1)). (3.10) 

N z (r\ — r 2 ) 1 

The left-hand side of (3.1) is usually called the wide (or smoothed) version of the 
eigenvalue density correlation function and the expression in the right-hand side 
of (3.10) is derived by various methods for different random matrix ensembles 

©BIS- 
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In Section 6 we study S(zi,Z2) with the spectral parameters zi,z% given 
by (3.8). Now the singularity of Q(zi,Z2) is determined by convergence of 
1 — vwiW2Uf(p) to zero. This convergence depends on the behavior of uf(j>) 
around the origin p = 0; that is why the rate of decay of u(t) at infinity dictates 
the form of the limiting expression for S in the local scale. 



3.2 The method and short proof of semicircle law 

We prove Theorem 3.1 in Sections 4 and 5. We are based on the moment 
relations approach for resolvents of random matrices proposed and developed 
in |2l], ^2[ ^8|. This technique is proved to be rather general, powerful and 
applicable to various random matrix ensembles. We use a modified version of 
this approach needed to study rather complex case of band random matrices. To 
make the subsequent exposition more transparent, let us describe the principal 
points of this method in application to the simplest case represented by GOE 
(2.6). 



3.2.1 Families of averaged moments 

In the early 70s F.Berezin observed Q that regarding correlation functions of 
the formal density of states Pn{8) = cr' N {0) 

Pi N \Q k ) = E PN (9 1 )--- PN (d k ), 

Bfc = (#1, . . . , 9k), one can derive for them a system of relations that resembles 

equalities for correlation functions of statistical mechanics. In this system P^ 

is expressed via sum of Pj^l, Pfc+i an d some terms that vanish in the limit 
N — > oo. This can be rewritten in the vector form 

=Po + BP^ +ft N \ 

with certain operator B and vector tf> such that ||S|| < 1 and || = o(l) in 

appropriate Banach space. These properties prove existence of lim jy— »ooP^ — 
P; the special form of B implies that the limiting P is nonrandom with the 
components Y\p{0j). 

This approach has got its rigorous formulation on the base of the resolvent 
approach used first in the random matrix theory in the pioneering work p4f . 
Regarding the resolvent Gn, the main subject is goven by the infinite system 
of moments 

k 

L[ N \Xk,Y k ;Z k )=El[G N (x j ,y j ;z j ), keN. (3.11) 
4=1 

The technique proposed in pl| , p8| and developed in p2| has been employed in 
the study of eigenvalue distribution of various ensembles of random operators 



and random matrices pfl, 26 
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In present paper we use the moment relations approach in its modified ver- 
sion. The main observation here is that often it is sufficient to study asymptotic 
behavior of and instead of the whole infinite family of the moments 
(3.11). This considerably reduces amount of computations and makes the proofs 
more transparent. To explain the principal steps of the proofs of Theorems 3.1 
and 3.2, let us present here the short proof of the semicircle law for GOE. 



3.2.2 Derivation of system of relations 

The main ingredients in the derivation of moment relations are the resolvent 
identity 

G'{z) - G{z) = -G(z) (H 1 - H) G'(z), (3.12) 
where G'(z) = {H 1 - z)- 1 , G{z) = (H - z)- 1 and equality 

E 7 F( 7 ) = E 7 2 E F'( 7 ), (3.13) 

where 7 is the Gaussian random variable with zero mathematical expectation 
and F(t), t G R is a nonrandom function such that all integrals in (3.13) exist. 
Equality (3.13) is a simple consequence of the integration by parts formula. 
Let us consider (3.12) with H' — An (2.6) and H = 0. We obtain relation 

N 

G N {x, y) = (5 xy G ^( x ' s ) A n(s, y), C = (3-14) 

8=1 

Regarding the normalized trace 

f N (z) = N-^GNfax) = (G N ) 
and using (3.13), we obtain relation 

N 



x,s=l 



One can easily find the partial derivatives with the help of (3.12). Remembering 
that H are real symmetric matrices, we have 

= - [G(x, s)G(t, y) + G(x, t)G(s, y)} (1 + S^y 1 . (3.16) 

Substituting (3.16) into (3.15), we obtain the first main relation for l[ N ^ 

E/iv = C + (v{VfN} 2 + ^ + 4 N \ (3-17) 

where 

= (vN-iE (G%) , and ^{ N) - (vE{f° N f N } 
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and we denoted by £° the centered random variable 



Clearly, (G}° = (G°) (here and till the end of the subsection we omit the 
subscript N in Gat). If one can show that two last terms in (3.17) vanish as 
N — > oo, then convergence E/at(z) — > w{z) will be proved. 

We estimate the term <pi with the help of two elementary inequalities that 
hold for the resolvent of a real symmetric matrix: 

\.f N (z)\ < \\G(z)\\ < \Imz\~ 1 

and 

||G 2 (z)|| < llmzp 2 . 
The last estimate implies that 

J2 \G(x, s)\ 2 = \\Ge x \\ 2 < \lmz\- 2 . (3.18) 

s 

Inequality (3.18) means that if z € A v , then (f>[ N ^ < vq^^N" 1 . 
3.2.3 Selfaveraging property 

To show that lim w ^oo V'i = Oj we prove that the variance of /jv vanishes 

Var/Ar = E|/^| 2 =0(iV- 2 ). (3.19) 

It is clear that 

VarjW = Ef° N f° N = Ef° N f N , 

where we denoted Jm = /n(z)- Applying (3.14) to the last factor /jv, we see 
that 

C N 

Vf° N f N =-± £ V{f° N G(x,s)A N (s,x)}. 

s,t=l 

The using (3.13) and (3.16), we derive relation 

E/^/at = Cv Vf° N fNfN + A N) + 4 N \ (3-20) 
where <p[ N) = (vN^Ef^ (G 2 ) and 

ip { 2 N) = 2(vN~ 2 ~E (G 2 G) . 
The useful observation is that 

E/^t/jv/jv = E/^/jvE/jv + E/^/^/jv- (3-21) 
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Using this identity and taking into account estimates (3.19), we derive from 
(3.18) that 

E/^/at < wry" 1 {Ef° N f N • E \f N \ + Ef° N f° N ■ \f N \} + 

v^N^E, \f^\ + 2vr]- 4 N- 2 . 
Taking into account that \f^\ 2 = /at/at ) we finally obtain 

E|/^| 2 =E/^/ w < 

1 /2 

2^- 2 E|/^| 2 +^ 3 ^ 1 (e|/ a V| 2 ) +2^- 4 iV- 2 . (3.22) 

This immediately implies (3.19) provided z£ A, (3.1). Obviously, t/>£ admits 
the same estimate. 

3.2.4 The semicircle law and further corrections 

Returning back to (3.17) and gathering estimates for (j)^ and ip{ N ^ , one can 
easily derive that if z £ A v (3.10, then limAr^oo gjv(z) = w(z) , with w(z) given 
by (2.10). Convergence of the Stieltjes transforms implies convergence of the 
corresponding measures. Thus the semicircle law is proved. 

It should be noted that relation (3.21) can be transformed into 

E/aV/at/at = 2E,f^f N Ef N + E/aV/aV/aV. 

Substituting this into (3.18), we see that 



Using the resolvent identity 



(3.23) 



G(z 1 )G(z 2 ) = - G{zi) G{Z2 \ G(zi) = (H - z^- 1 (3.24) 
zi - z 2 

and convergence of E/at(z), one can easily find the limiting expression for 
E(G 2 G). If one assumes that two last terms in (3.23) are values of the order 



o(N 2 ), then one arrives at (3.7) (see e.g. |20 for more details) 



4 Correlation Function of the Resolvent 

Our approach is to apply systematically the scheme of subsection 3.2.2 to get 
the leading term of the correlation function C^ n ' b '(zi, z%) (3.4). This term is 
expressed via the limit of the limE(G^ n,b ^(z)) = w(z) but we have to prove the 
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pointwise version of this convergence given by Theorem 3.1. This and other 
auxiliary propositions are addressed in subsections 4.1 and 4.2. In subsection 
4.3 we give the proof of Theorem 3.2 on the base of these statements. In what 
follows, we omit super- and subscripts n and b and do not indicate the limits of 
summation when no confusion can arise. 

4.1 Proof of Theorem 3.1 

Using relations (3.12)-(3.14) with obvious changes and repeating computations 
of subsection 3.2.2, we obtain relation 

EG(x,x) = C + CvVG{x,x)U G {x) +(v E [G(x, s)] 2 U(s, x), (4.1) 

\s\<n 

where 

U G (x) = J2 G(s,s)U(s,x) = ± J2 G(s,s)u( S -^). 

|s|<n M<™ 

Let us denote the average EG(x,x) by g(x) and rewrite (4.1) in the following 
form 

g(x) =( + (v 2 g(x) U g (x) + ^(x) + *(x), (4.2) 
where we denoted (cf. (3.17)) 

$W,(^E[G(v)] 2 ^) (4.3) 

\s\<n 

and 

tf(a;) = (vEG°(x,x)U G (x). (4.4) 
Let us consider the solution {r(x), \x\ < n} of equation 

r(x) =( + (vr(x)U r (x), \x\ < n. (4.5) 

Given z E A v (3.1), one can prove that the system of equations (4.5) is uniquely 
solvable in the set of A^-dimensional vectors {r} such that 

\\r\\ 1 = sup \r(x)\ < 2r]-\ rj = \Imz\ (4.6) 

\x\<n 

(see Lemma 4.1 at the end of this section). Certainly, r(x) depends on particular 
values of z, n and b, so in fact we use denotation r(x) — r n ^(x; z). 
The following statements concern the differences 

D n , b (x; z) = g n , b (x; z) - r nfi (x\ z), d nfi {x; z) = r„ ib (x; z) - w(z), 

where w(z) is given as a solution of (2.9). 



12 



Proposition 4.1. 

Given e > , there exists a number L = L(e) such that for all sufficiently 
large b and n satisfying (2.13) 

sup \d n .b(x; z)\ < e, ze A„ (4.7) 

xeB L 

with Bl given by (3.2). 
Proposition 4.2. 

If z e A v (3.1) and (2.13) holds, then 

sup \D n ,b(x; z)\ = o(l), n, 6 ^ oo. (4.8) 

\x\<n 



Theorem 3.1 follows from (4.7) and (4.8). Under the same conditions one 
can find V > L such that 

s'Kk-" 11 ^' (4 - 9) 

Relation (4.9) follows from (3.3) added by (4.6), a priori estimate 

sup \g(x)\ < (4.10) 

|x|<n \ lmz \ 

and observation that V has to satisfy condition u(L — L') < e. 
Proof of Proposition 4-1 

Let us consider the constant function w x (z) = w(z) satisfying (2.10) that we 
rewrite in the following form similar to (4.5) 

w x (z) = ( + (vw x (z)- ^2 bS xt w t (z), \x\ < n. 

\t\<n 

Subtracting this equality from (4.5), we derive that d(x) = d n ,b(x;z) verifies 
equality 

d(x) = (vd(x)U r (x) + (vw(z)U d (x) + (vw 2 (z) [P b + T{x)\ , 

where 

p ^- h j2 u (i)- r <s)d 8 (4.H) 

tez v ' J -°° 

and 

't 



Tn , b ( X ) = lj2 u (^-)-- b Y, u ( 



\t\<n x / tez 



bl- ^ 
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It is clear that \Pb\ = o(l) as b — > oo. Indeed, one can determine an even 
step- like function Md(i), t e R such that 

u d {t) = ^<\)\h^lh){t),t>Q. 

Then Ud(t) < u(t) and Ud(t) — > u(t) as — > oo and the Beppo-Levy theorem 
implies convergence of the corresponding integrals of (4.10). 
Taking into account equality 

r(x) 



l-v(U r (x) ' 

we can write that 

d(x) = vwr(x)U d (x) + vw 2 r(x) [P b + T rhb (x)} , 

where we denoted w = w{z). This relation, together with estimates (4.6) and 
\w(z)\ < |Imz| _1 , implies inequality 



sup \d{x)\ < r sup \d{x)\ + sup |T„ ;b (a;)| + P b < 

xeB L \xEBl-! xeB L J 

tJ^t 3 ( sup |T„ i6 (a;)|+n J + 

t l sup \d(x)\ , (4-13) 



3=0 



iEBi-j / \x\<n 



where r < vr] 2 < 1. It is clear that due to monotonicity of u(t), one gets 

2 00 /f\ r°° 
sup |T„. 6 (a;)| < sup |T„, 6 (a;)| <- V u - < 2 / u(s)ds. 
xeB L+1 xeB L o t J^Z Lb \oj J L 

Given e , one can find such a number k that r fc < e/4. Than we derive from 
(4.13) that 

k 

sup \d(x)\ < t^V su P |r„,b(x)|+2rP b + e/4< 

x£B L ._ Q xeB L _j 

r(fc + l) sup |T n , 6 (;r)| + 2rP b + e/4. 

Now it is clear that (4.9) holds for sufficiently large L and b. 
Proposition 4.1 is proved. 

Proof of Proposition 1^.2. 
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Subtracting (4.5) from (4.2), we obtain relation for D(x) = D n fi{x) 

1 



D{x) = (vD(x)U r (x) + (vg(x)U D (x) + (v 
that can be rewritten in the form 



-$(x) + *(x) 



D(x) = vg(x)r(x)UD{x) + vr(x) 



i$(x) + #(x) 



Regarding this relation as the coordinate form of a vector equality, one can write 
that 



D = v(l -W {9 ^ 



b 



where we denote by W^ 9 '^ a linear operator acting on vectors e with components 

e(x) as 

(x) = vg(x)r(x) ^ e(s)U(s,x) 

\s\<n 



and vectors 



= r{x)4> n , b {x), ¥ r \x) = r(a;)*(a;). 



It is easy to see that if z e A,,, then the estimates (4.6) and (4.10) imply 
inequality 



<- 2 < 1/2. 



Thus, to prove Proposition 4.2, it is sufficient to show that 



sup 



5>[G(x, s )] 2 M ( 



s — X, 



Oil), zeA v 



and 



sup'E\G ix 1 x)Uaix,x)\ =o(l), ze A v . 



(4.14) 



(4.15) 



(4.16) 



Relation (4.15) is a consequence of the bound (2.2) and inequality (3.18). Re- 
lation (4.16) reflects the selfaveraging property of G^ n ' b \ This question is ad- 
dressed in the next subsection. It should be noted that (4.16) will be proved 
independently from computations of this subsection. Assuming that this is done, 
we can say that Theorem 3.1 is proved. We complete this subsection with the 
proof of the following auxiliary statement. 

Lemma 4.1. 

Equation (4-5) has a unique solution in the class of vectors satisfying condi- 
tion (4-6). 



Proof. 
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Let us consider the sequence of TV-dimensional vectors {f^ k \ k e N} deter- 
mined by relations for their components 

r (fe+1) (x) = C + (vr (k) {x)U rW (x), r (1) (x) = £ M < n. 

Then it is easy to derive that if r^ fe ' satisfies (4.6) and z G A n (3.1), then 
f^ k+1 ^ also satisfies (4.6). The difference Xk+\{x) = r^ k+1 \x) — A k \x) satisfies 
relations 

Xk+i{x) = (v Xk {x)U rW (x) + Cur (fe - 1} (x)U Xk (x). 
Obviously, < a Hxfclli with a < 1 provided zeA,. Lemma is proved. □ 

4.2 The variance and selfaveraging property 

The asymptotic relation (4.15) is a consequence of the fact that the variance of 
G(x, x) 

Var(G ( "< 6) ) = E |(G (n < b) )° 2 = E {{G {n ' b \z))° {G {n ' b) {z))°} 

vanishes as n, b — > oo. Instead of the direct proof of (4.15), we prefer to present 
the whole list of more general statements needed in studies of the correlation 
function of G. All of them can be proved independently of the Theorem 3.1 
without use of its statement. 

We start the list with the following three relations that concern the moments 
of diagonal elements of G. 

Proposition 4.3. 

If z e A,, (3.1), then the estimates 

sup E|G°(a;,a;;z)| 2 = O^ 1 ), (4.17) 

|x|<n 

sup V\U° G (x)\ 2 -0(6- 2 ), (4.18) 

\x\<n 

and 

sup E\U^(x)\ 4 = 0(b- 4 ), (4.19) 

\x\<n 

hold. 

The following statement concerns the mixed moments of variables G° (x, x; z) 
and their sums. 

Proposition 4.4. 
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If z € A r) , then relations 

sup \EG°(x,x)U° G (y)\=0(b- 2 ), (4.20) 

|x|,|i/|<n 

sup \-E(G }G°(x,x)\=o(n- 1 b- 1 +b- 1 [Var(G)} 1/2 ) , (4.21) 

and 

sup |E(G )G°(x,a;)^(y)| =0(n- 1 fe- 2 + 6- 2 [Var(G)] 1/2 ) (4.22) 

\x\,\y\<n V ' 

are true in the limit n,b — > oo. 
Finally, we formulate 



Proposition 4.5. 

//zeA,, i/ien relation 



sup 

|x|<n 



E (G^^s)] 2 ^,*) 



O^-^ + ft" 2 [Var(G)] 1/2 ) 

(4.23) 



is true in the limit n,b — > oo. 



Let us not that the estimates (4.21)-(4.23) admit also the estimates in terms 
of n and b that do not involve the variance of (G). However, derivation of the 
estimates would take more place and taime and we restrict ourselves with the 
forms presented. It will be shown later that Var (G) = 0(n~ 1 b~ 1 ). This fact 
together with the restriction (2.12) implies for (4.22) and (4.23) that 

1 1 1 



b 2 ^/nb nb 

that is sufficient for us. We prove Propositions 4.3-4.5 in Section 5. 

4.3 Toward the correlation function 

Let us have a more close look at the correlation function 

C n , b (z u z 2 )=v{(G( n > b \z 1 )) (G( n > b Hz 2 ))°} 

We follow the scheme described at the end of subsection 3.2 and introduce 
variables Gj(x,y) = G^™'^ (x, y; Zj), j — 1,2. To study the average 

E(G 1 )G 2 (a ; ,a ; ) = J R 12 (a ; ), 
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we apply to G2{x,x) the resolvent identity (3.12) and obtain relation 



Rx2{x) = -& J2 V{(G 1 )G 2 (x, S )a(s,x)}^/U(s,x), 

\s\<n 

where (2 = — z 2 . We compute the last mathematical expectation with the help 
of formulas (3.13) and (3.16) and obtain equality (cf. (4.1)) 

Ri2(x) = (2vRi2(x)U g2 (x) + (2vlI Rl2 {x)g 2 (x) + 
2C2VN- 1 EG?(x, s)G 2 (x, s)U(s, x) + ( 2 v [Q 12 (x) + T 12 (x)} , 

s 

where we denoted g 2 (x) = ~EG(x, x;z 2 ), 

Ug 2 {x) = 92{S)U(S,X), 
\s\<n 

Ur 12 (x) = J2 Ri2(s)U(s,x), 

\s\<n 

e 12 (x) = e J (Gi) J2 s )] 2 u (*> x ) I . 
I l«l<„ I 



and 



r 12 {x)=E{{Gl)U° G2 {x)G° 2 {x)}. 

Using denotation 



«■<•) - wfcw (4 - 24) 

we obtain the following relation for R 12 

Ri2(x) = vq 2 (x)g 2 (x)U Rl2 (x) + ]T F 12 (x, s)U(s, x)+ 

S 

vq 2 {x)[e 12 {x) + T 12 {x)l (4.25) 

where we denoted 

F 12 (x,s) = EGt(x,s)G 2 (x,s). 

The terms 6 and T can be estimated with the help of Propositions 4.3-4.5. 
As we shall see in the next subsection, they do not contribute to the leading 
term of Ri 2 . To obtain the explicit expression for the this leading term, it is 
necessary to study in detail the variable F12 . Now let us formulate corresponding 
statement and the auxiliary relations needed. 
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Proposition 4.6. 

If z e A v , (3.1), then for arbitrary positive e and large enough values of b 
and n (2.13) there exists the set Bl (3.2) with L such that 

sup 

x<£B L 



The proof of Proposition 4.6 is based on the similar statement formulated 
for the product G\G 2 . 

Proposition 4.7. 

Given positive e, there exists such L that relations 

sup 

xeB L 

(4.27) 

and 

sup 

x£B L 

hold for all k > 1, all Zi G A v and large enough values ofb. 

Remark. In the case when z\ ^ Z2, relation (4.28) can be derived from the 
resolvent identity (3.24) with the help of the convergence (3.3) and the explicit 
form of w{z) (2.9). 

We prove Proposition 4.6 in the next subsection. Relations (4.27) and (4.28) 
will be proved in Section 5. 

4.4 Proof of Proposition 4.6 and Theorem 3.2 

Let us assume that relations (4.27) and (4.28) are true and show that under 
conditions of Theorem 3.2 the leading term of i?i 2 is of the order 0(n~ 1 b~ 1 ) 
and terms 0i 2 and T 12 of (5.2) do not contribute to it. We rewrite (4.25) in 
the form 

R 12 (x) = vg 2 (x)q 2 (x)U Rl2 (x) + 2vq 2 (x)N- 1 [F 12 U] (x,x)+ 

vq 2 (x)[G 12 (x) + r 12 (x)}. (4.29) 
Let us denote r\ 2 = sup^^^ |i?i 2 (x)|. 
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b[F 12 U](x,x) 



w\w 2 



u F (p) 



2tt 1 - vw 1 Jr [1 - vwiw 2 u F (p)] 



T dp 



< e. 



(4.26) 



2ir J R 1 - vw 1 w 2 u F (p) 

\s\<n 



EG 1 {x,s)G 2 (x,s) - - 

|s|<n 



WlW 2 



VW\W 2 



< e 



(4.28) 



Taking into account U(x,y) < u/b (2.2) and using inequalities of the form 
(3.19), it is easy to see that if z% e A ?) , then 

(\ 1/2 / v 1/2 

Regarding this estimate and relations (4.22), (4.23) we easily derive from (4.29) 
inequality (cf. (3.22)) 

^ v CI, 

t-12 < ~^r 12 + — + -nrVri2 
t] z bn b A 

with some constant C. Since r\ 2 is bounded for all z£ A„ then 

Now condition (2.12) implies that r\ 2 = 0(l/nb) and therefore the general form 
of (3.5) is demonstrated. 

Substituting (3.5) into the estimates (4.22) and (4.23), we obtain that 

1161211,= and 11^=0(1). 

Thus, these terms of (4.29) do not contribute to the leading term of R 12 . To 
find this term in explicit form, we need the result of Proposition 4.6. 

Proof of Proposition 4- 6. 

Regarding F\ 2 {x, y) = E G\ (x, y)G2(x, y), we apply to G 2 the resolvent iden- 
tity (3.12). Computations similar to those of subsection 3.2.2 lead us to equality 

F 12 (x,y) = <; 2 5 xy EGj(x,x)+( 2 v[t 12 U] (x,y) EG 2 1 (y,y)+ 

C2V {[F 12 U] (x,y) 9l (y) + F 12 (x,y)U [g2] (y) + T(x,y)} , (4.30) 

where 

t 12 (x,y) = ET 12 (x,y) = EG 1 (x,y)G 2 (x,y), 
and the vanishing terms are denoted by T = Ti + T 2 + T^: 

r x (x, y) = ^E {G^x, y)G\{s, y)G 2 {x, s) + 2G\{x, y)G 2 (s, y)G 2 (x, a)} U(s, y), 

s 

T 2 (x, y) = E { [T 12 U] (x, y) [G\{y, y)]°) + E{[F 12 U] (x, y)G° 1 (y, y)} , 

and 

r 3 (x,y)=EF 12 (x,y)U( G2] (y). 
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Indeed, it is easy to show that 

sup \Tj(x,y)\ = 0(& _1 ), z 1 ,z 2 eA r] . 



(4.31) 



This can be done with the help of inequality (3.19) and relations (4.17), (4.18), 
and (4.23). 

Using definition of 52(2;) (4.24), we rewrite (4.30) as 



where we denoted 



Fi2(x,y) = vg 1 (x)q 2 (y) [F 12 U] (x,y) + 
R {1) {x,y) + R {2 \x,y)+vf{x,y), 

R {1 \x,y) =q 2 {x)EGt(x,x)S xy , 
R^(x, y) = vq 2 (y) [t 12 U] (x, y) EGj(y, y) 



(4.32) 

(4.33a) 
(4.336) 



and f(x,y) = T(x,y)q 2 (y). Let us note that |i?(l)| < r^ 3 and |i?( 2 )| < vr}~ 5 for 
z l e A v . 

Let us determine the linear operator W that acts on N x N matrices F 
according to the formula 



[WF}(x,y)=v 9l (x) 



£ F(x,s)U(s,y) 

\s\<n 



Q2 (y)- 



The a priori estimates |<7i(a;)| < |Imzi| 1 , and |?2(a;)| < jlm^l 1 imply in- 
equality (cf. (4.14)) 



v 1 



(4.34) 



||W||(i,i) < -J < o. Zi6A„ 

where the norm oi N x N matrix A is determined as ||^4||(i ,i) = sup x y \A(x,y)\. 
This estimate is verified by direct computation of 1 1 T4^^4_| j (1,1) with || ^4|| (1.1) = 1- 
Then (4.32) can be rewritten as 

00 

F 12 (x,y) =vJ2 [ wm ( R{1) + r(2) +vt)] (x,y). (4.35) 

m=0 

The next steps of the proof of (4.26) are very elementary. We consider the 
first M terms of the infinite series and use the decay of the matrix elements 
U(x,y) = U^ bS) {x,y) . Indeed, if one considers (4.33) with x and y taken far 
enough from the cndpoints —n, n, then the variables gi(s), q 2 (t) enters into 
the finite series with s and t also far from the endpoints. Then one can use 
relations (3.3) and (4.9) and replace gi and q 2 by the constant values wi and 
w 2 , respectively. This substitution leads simplifies expressions with the error 
terms that vanish as n, b — > 00. The second step is similar. It is to show that we 
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can use Proposition 4.7 and replace the terms and R^ of the finite series 
of (4.33) by corresponding expressions given by formulas (4.27) and (4.28). 

Let us start to perform this program. Taking into account the estimate of T 
and using boundedness of terms R^ 1 ' and R( 2 ', we can deduce from (4.35) that 



M 

bJ2 F Mx,s)U(s,x) =bvJ2 [\V m (R {1) +R {2) )U\ (x,x) + A^(x,x), (4.36) 



m=0 



where M is such that given e > 0, \A^(x, x)\ < e for large enough b and n. 
Now let us find such h that the following holds 

u(t) < e,y\t\ > h, and / u(t)dt < e. 

J\t\>h 



We determine matrix 



K ' yj 1 0, if \x - y\ > bh 



and denote by W corresponding linear operator 



[WF}{x,y)=v 9l {x) 



£ F(x,s)U(s,y) 

\s\<n 



Q2{y)- 



Certainly, W admits the same estimate as (4.34). 

Given e > 0, let L the largest number among those required by conditions 
of Propositions 4.1 and 4.7. Let us denote by Q the first natural greater than 
(M + k)h. Then one can write that 



bv 



M 

[w m {R {1) +i? (2) )t/] (x,x) = 



m=0 



where 



M 

bvJ2 [{vW) m {R (1) +i? (2) )L>] {x,x) + A^(x,x) 



m=0 



sup \A^ 2 \x,x)\ < e, as n, in oo. 

x£B l+ q 



(4.37) 



The proof of (4.37) uses elementary computations. Indeed, A^(x,x) is repre- 
sented as the sum of M + 1 terms of the form 



bv m+1 Y [gi(x)] m F(x,x 1 )U(x,xi)q 2 (xi) ■ ■ ■ U{x m -2,x m - 1 )q 2 {x m -i)x 

\xi\<n 
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+ i? (2) ](x m _i, x m )U(x m , x), 

where the sum is taken over the values of Xi such that \xj — Xj+i\ > bh at least 
for one of the numbers j < m. 

Now remembering the a priori bounds for i?W and R^ 2 \ estimates like (4.13) 
and taking into account the diagonal form of one obtains the following 

estimate of A^ 2 ^ by two terms 



M * 

sup |A (2) (a;,a;)| < ^ ^ 2m+3 U(x,Xi)U(x 1 ,x 2 ) ■ ■ ■ U(x m ,x)+ 

\xi\<n 



\x\<n 



m=0 



M v m+l 



E 



2m+5 



^ U(x, Xi) - ■ ■ U(x m -2,x m -i)U(x m ,x). 



(4.38) 



m=0 



\xi\<n 



Regarding the first term in the right-hand side of (4.38) and assuming that 
\xj — xj + i\ > bh, one can observe that for large enough b, n 



Indeed, 



J2 U j (x, Xj )sU m - j (x j+1 ,x) < e. 

\xj\<n 

U(x,x 1 )U(xi,x 2 ) ■ ■ ■ U(xj-i,Xj) < 

\xi\<n 



_ r poo 

U{x,xi)U{xi,X2) ■ ■ ■ U(xj-i,xj) < / u(t)dt + 

„ r*7 \_J — OO 



XiEZ 



u(0) 



< {i + u/by. 



Let us also mention here that given e > 0, one has for large enough n, b that 



sup 

x£B l+ q 



\s\<n 



U 3 {x,s) - 1 



(4.39) 



where j < M. This follows from elementary computations related with the 
differences (4.11) and (4.12) that vanish in the limit l«k<«. 

Similar but a little more modified reasoning can be used to estimate the 
second term in the right-hand side of (4.38). Now one can write that 



M 



sup |A (2) (x,a;)| < 2e ^ 



rnv 



m+l 



|x|<n 



„2m+2 



< e. 



m=0 



Regarding the right-hand side of (4.37) with x G B l+ q, one observes that 
the summations run over such values of Xi that |ar — a?i j < bh, \xi — x i+ i \ < bh, 
and thus xj £ Bl for all j < k + m — 1. This means that we can apply relations 
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(3.3) and (4.9) to the right-hand side of (4.37) and replace gi(x) by w(zi), q2(x) 
by w(z 2 ). We derive from (4.36) that 

M 

(F 12 U k ){x,x) = bviD 2 ^(vw 1 w 2 ) m {U m )(x,s) \rW + R^A (s,t)U(t,x)+A^ (x,x) 
with 

sup |A (3) (x,a;)| < As. 

xeB L+ Q 

Finally, applying Proposition (4.7) to the expressions involved in R and taking 
into account that 

sup \bU m+1 (x,x) - A- / u% +1 (p)dp\ < e, (4.40) 
xeB L+Q ^ J 

we obtain equality 

2 M r 

(F 12 U)(x,x) = l-^_^ vw ^) m / «F +1 b)dp+ 



m=0 

fr^E^r / 1 fi?+1(p) - d P+ A^) ( ,,,) (4.4i) 

Z7T 1 — utyf ^— / 1 — vw\w 2 uf 

with 

sup |A^ 5 ^(x,a;)| < £ b,n^oo. 

xeB L+Q 

Passing back in (4.41) to the infinite series and simplifying them, we arrive 
at the expression standing in the right-had side of (4.26). Proposition is proved. 
□ 



Let us complete the proof of Theorem 3.2. Remembering estimate (4.14), 
we can iterate relation (4.29) and obtain that 

Ri2(x) = ^^f^[(W^rf 12 ](x)+o(l/nb), 

where we denoted f\ 2 (x) = bq 2 {x)[F 12 U]{x , x). Regarding the trace 
1 £ R u {x) = ± Y, ^12(^(1 + 0(1)), 

ja:j<n x£Bl 

and repeating the arguments of the proof of Proposition 4.6 presented above, 
we can write that 

M 

Ri2(x) = Yl E( 6F i2t/)(t,i)(^ 2 2 t/r(i,x) + a^(x,x) 

m=0 t 
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with sup xeBi |A( 6 )(x,a;)| < vep' provided n, b — > oo (2.12). Finally, observing 
that {bF 12 U){t,t) asymptotically does not depend on t (4.26), we arrive, with 
the help of (4.39), at the expression (3.6). Theorem 3.2 is proved. 



5 Proof of auxiliary statements 

Proof of Proposition 4-3 

Let us consider the average E G°(x, x)G 2 (y, y) and derive for it, with the 
help of formulas (3.12), (3.13) and (3.16) equality 

E G\(x, x)G 2 (y, y) = C 2 «E G\ (x, x)G 2 (y, y)U G2 (y)+ 
( 2 vJ2VG°(x,x) [G 2 (y }S )} 2 U( S ,y) + 

S 

2C2vY,VG 1 (x,s)G 1 (y,x)G 2 (y,s)U(s,y). 

s 

Applying to the first term of this equality the analogue of identity (3.21) 
and using q 2 (x) (4.24), we obtain that 

EG° 1 (x,x)G 2 (y,y) = vq 2 (y)EG° 1 (x,x)G 2 (y,y)U G2 (y)+ 
vq 2 (y) J2 EG °i (*> x ) I G 2 (y, s)f U(s,y)+ 

S 

2vq 2 (y) E Gi(x, a)Gi(y, x)G 2 (y, s)U(s, y). (5.1) 

s 

We multiply both sides of this relation by U(x,t) and sum it over x; then we 
get 

E U° Gi (t, t)G 2 (y, y) = vq 2 (y)E U° Gi (t)G 2 (y, y) U° G2 (y)+ 
vq 2 (y) E U Gl (t) [G 2 (y, s)} 2 U(s, y)+ 

S 

2vq 2 (y) E G^x, s)G 1 (y, x)G 2 (y, s)U(s, y)U(x, t). (5.2) 

s 

Regarding G\{y,-)U{-,t) and G 2 (y , -)U (• , y) in the last term as vectors in 
iV-dimcnsional space, we derive from estimate (3.19) that 



J2 E Gx(x, s)G 1 (y, x)G 2 {y, s)U{s, y)U{x, t) 



< 



1/2 



1/2 



\GA\[Y,\Gx{v,x)U{x,t)\ 2 ) [^\G 2 (y,s)U(s,y)f 



(5.3) 
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Inequality (3.18) implies that the right-hand side of (5.3) is bounded by b~ 2 n~ 3 . 

Let us multiply both sides of (5.2) by U(y,r) and sum them over y. Then 
one obtains a relation that together with (3.18) and (5.3) implies the following 

estimate for variable M\ 2 = sup x f E | U Gi (x) | ) : 

M 12 < vrf 2 M 12 + vq- 3 b- 1 x /M^ + 2vr]~ 4 b~ 2 . 

This proves (4.18). 

Now (4.17) follows from (4.18) and relation (5.1). 
To derive estimate (4.19), let us consider the variable 

E U° Gl ( Xl )U G2 (x 2 )U G3 (x 3 )U° Gi (x 4 ) = E [U° Gl (x)U° G2 (x)U° G3 (x 3 )]° U Gi {x 4 ). 

Let is denote T = U G U G U G3 and and M(x t , x 2 , x 3 , t) = F,T°G 4 {t,t). We 
apply to G 4 (t,t) resolvent identity (3.14) and obtain relation 

ET°G 4 (M) = v( A ET°G i (t,t)U Gi {t) + 
v( 4 ET°Y,[G4(s,t)} 2 U( S ,t) + 

S 

< 4 Yl ^U Gi (x i )U G .(x j )J2G k (y,s)G k (t,y)U(y,x k )G 4 (t,s)U(s,t). (5.4) 

{i,j,k) x,s,t 

Repeating previous computations and applying similar estimates, we obtain 
inequality 

\Y,M(x u x 2 ,x 3 ,t)U(t,x 4 )\ < -E\TU Gi (x 4 )\ + -E\T\E\U Gi (x 4 )\ + 
t ^ ^ 

^ b nT\ + ^nu G M)u° G Mi)\- ( 5 - 5 ) 

Here we have applied inequalities (3.18) and (5.3) to the last two terms of 
relation (5.4). Now it is clear that (5.5) implies (4.19). Proposition 4.3 is 
proved. □ 

Proof of Proposition 4-4- 

Estimate (4.20) follows from relation (5.2) and estimate (4.18). Regarding 
(5.1) and summing it over x, one can easily derive (4.21) with the help of the 
arguments used to prove (4.18). 

Let us turn to the proof of (4.22). To do this, let us consider the variable 

K(x, y) = E (G°) G°(x, x)U° G {y) = E [(G°) U G (y)}° G(x, x) 
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and apply to the last expression resolvent identity (3.12) and formulas (3.13) 
and (3.16). We obtain equality that can be written in the following form with 
denotation R=(G°) U G (y) 

ER°G(x,x) =CvER°G(x,x)U G (x)+ ^ iu(x,y), (5.6) 

i=l,2,3 

where 

Kl (x, y) = (v E R°G(x, s)G(x, s)U(s, x), 

S 

k 2 (x, y) = 2(v ^2 E (G°) G(t, s)G(x, t)u 2 b {t, y)G(x, s)U(s, x), 

s,t 

and 

k 3 (x, y) = 2CVN- 1 J2 s)G(x, t)U° G {y)G(x, s)u 2 (s, x). 

s,t 

Let us use identity 

~E,R°XY = ~ERX°-E,Y + ERY°EX + ERX°Y° - ~E R~E X°Y° . 
and can rewrite (5.6) in the form 

E R ° G[X ' X) = l-lf{x)g{x) [E RU °o^ G ° ^ x ) - E < G °> U °g{v)VG (x)U° G (x)\ + 

1 - vg(a:)s(x) .J^ 3 

Taking into account relation (4.18), inequalities (3.18) and (5.3), we obtain 
that 

\Ki(x,y)\ < 2 V - 2 b- 2 (Var (G)) 1/2 fori = 1,2 

and 

\K 3 (x,y)\ <2r 1 - 3 b- 2 N- 1 . 
Using them, we derive from (5.7) inequality 

|if( 2 ;, 2/ )|<2r / - 1 (Var (G)) 1/2 | (e \U g (x)\ 4 ) ^ + b~ 2 (e \U g (x)\ 2 ) ^ 

2r 1 - 1 b- 2 (Var (G) ) 1/2 + 2r]- 2 b- 2 N~ x . 
This leads to estimate (4.22). Proposition 4.4 is proved. □ 

Proof of Proposition 4-5. 

This proof of the estimate (4.23) is the most cumbersome. Here we have to 
use the resolvent identity (3.12) twice. However, the computations are based 
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on the same inequalities as those of the proofs of Propositions 4.3 and 4.4. 
Therefore we just indicate the principal lines of the proof and do not present 
the derivations of estimates. 

To compute the mathematical expectation 

EM(x, S ) = E(G° 1 ){G 2 (x, S )}\ 

let us apply to G 2 (x,s) the resolvent identity (3.12). We obtain equality 

E M(x, s) = C 2 ^E (G°) G 2 (x, x)- 

C 2 E (G?> G ^ x ' s)G 2 (x, t)a(t, 8 )y/U%8). (5.8) 
t 

Relation (4.21) implies that the first term of the right-hand side of (5.8) is the 
value of the order indicated in (4.23). Let us consider the second term of (5.8). 
We compute mathematical expectation with the help of relations (3.13) and 
(3.16) and obtain expression 

5 

c 2 e (g°) y, G ^ s ) G i( x > *)«(*. s)VW^j = E e *( x ' s )' ( 5 - 9 ) 

t i=l 

where 

&i(x, s) = < 2 E (G?) G 2 (x, s)G 2 {x, s) EU G (s), 
e 2 (x, a) = v( 2 E (Gf ) G 2 (x, s)G 2 (x, 

9 3 (x, s ) = ^E^G?( S ,t)(7(t, S )G 2 (x, s )G 2 (x,t), 
t 

04(x, a) = < 2 E (Gf) ]T [G 2 (x, t)] 2 U(t, s)G 2 (s, s), 
t 

and 

Q 5 (x,s) =2< 2 E (£?J)5^G2(a;,s)G2(a;,t)G?2(s ) t)I7(t,a). 
t 

9i is of the form ^EMfi, s)E[/g(s) and can be put to the right-hand side of 
(5.9). The terms 2 and 83 are of the order indicated in the right-hand side of 
(4.23). This can be shown with the help of estimates of the form (5.3). 

Regarding 64, we apply the resolvent identity (3.12) to factor G 2 (s, s). Re- 
peating the usual computations based on (3.13) and (3.16), we obtain that 

G 4 (x, s) = vQl Y EM(x, t)U(t, s) + v( 2 4 (x, s)EU G2 (s) + Q{x, s), (5.10) 
t 

where f2 gathers the terms that are all of the order indicated in (4.23). This can 
be verified by direct computation with the use of estimates (4.18), (4.21), and 
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(4.22). Not to overload this paper, we do not write down the terms constituting 
SI and do not present their estimates as well. Relation (5.10) is of the from that 
leads to the estimates needed for ^2 S EM(x, s)U(s,x). 

Regarding Q 5 (x,s), we apply (3.12) to G2{s,t) and obtain, after the use of 
(3.13) and (3.16) that 

6 5 (x, s) = 2v(i^EM{x, s) + v(0 5 (x, s)EC/ G2 (s) + ft'(x, s), (5.11) 
o 

where Q'(x, s) consists of the terms that are of the order indicated in (4.23). 
The form of (5.11) is also such that, being substituted into (5.9) and then into 
(5.8), it leads to the estimates needed. This observations show that (4.23) is 
true. 

Proof of Proposition 4-7. 

We prove relation (4.27) with k = 1 because the general case does not 
differ from this one. To derive relations for the average value of variable 
ti 2 {x,y) = EGi(x,y)G 2 (x,y), we use identities (3.12)-(3.14) and repeat the 
proof of Proposition 4.6. Simple computations lead us to equality 

h 2 (x,y) = gi(x)( 2 S xy + ( 2 v 2 t 12 (x,y)U g2 {y)+ 

4 

( 2 v 2 ]T tia (x,s)U(s, y) 9l (y) + ( 2 v 2 ]T T j (x, y), (5.12) 

j=i 

where 

Tifo y) = Ej2Gi(x, y)G 2 {x, s)G 2 {y, s)U(s, x), 

S 

T 2 (x,y) = E^2G 1 (x,y)G 1 (s,y)G 2 {x,s)U(s,x), 

S 

T 3 (x,y) = EG 1 (x,y)G 2 (x,y)U^ 2 (y), 

and 

T 4 (x, y) = EG°(y, y) £ G 1 (x, s)G 2 (x, s)U(s, y). 

S 

It is easy to see that inequality (4.16) implies estimates 

sup |Ti(as, y)\ < 6" V 3 , sup | V T^x, y)\ < 6~ V 3 - 

x,y x 

>« y 

The same is valid for T 2 (x, y). Similar estimates for T^(x 1 y) and T^(x 1 y) follow 
from relations (4.17) and (4.18). Thus, (5.12) implies that 

h 2 {x, y) = gi(x)q 2 (x)6 xy + vg 1 (y)q 2 (y) [t 12 U] (x, y) + A(x, y), (5.13) 
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where 



sup |A(:r, y)\ = o(l) and sup 

x,y x 



o(l) (5.14) 



in the limit n, b — > oo (2.12). 

We rewrite relation (5.13) in the matrix form (cf. (4.35)) 

_ , oo 

h2 = (i - W^) [Diag( 5l(?2 ) + A] = [W^] m (Diag( 5l(72 ) + A) 



m=0 



(5.15) 

Now we can apply to (5.15) the same arguments as to (4.35). Replacing g\{x) 
and q 2 (x) by Wi and W2, respectively, we derive from (5.14) that for x € Bl + q 

M 

tv»{x, s)=J2 (wiw 2 ) m+1 [U m ] (x, s) + o(l), n, b -> 00. (5.16) 

Multiplying both sides of (5.16) by U(s,x) and summing it over s, we obtain 
relation 

M 

t 12 {x,s)U(s,x) = (wiw 2 ) m+1 [U m+1 ] (x,x)+o(l), n,6- 00. (5.17) 

\s\<n m=0 

Now convergence (4.40) implies relation that leads, with M replaced by 00, to 
(4.27). 

To prove (4.28), let us sum (5.16) over s. The second part of (5.14) tells 
us that the terms A remains small when summed over s. Thus we can write 
relations 

M 

ti 2 {x, s)=Y^ {w lW2 ) m+l i um ] 0> s ) + n,b->oo. (5.18) 

s m— s 

Taking into account estimates for terms (4.11) and (4.12), it is easy to observe 
that convergence (4.39) together with (5.18) implies (4.28). □ 

6 Asymptotic properties of S(zi,Z2) 

In the last decade, the main focus of the spectral theory of random matrices is 
related with the universality conjecture of local spectral statistics put forward 
first by F. Dyson [ p^[ . This problem is addressed in a large number of papers 
where various random matrix ensembles are studied using different approaches 
(see e.g. the review [|l6|). The best understood are the Gaussian Unitary En- 
semble (GUE) and its real symmetric analogue GOE (see (2.5)). The proba- 
bility distribution of these ensembles are invariant with respect to the unitary 
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(orthogonal) transformations. This leads to the fact that the joint probability 
distribution of eigenvalues of these ensembles does not depend on the distribu- 
tion of eigenvectors and is given in explicit form p5| . This allows one to use 
the powerful technique of the orthogonal polynomials that provides a detailed 
information of the spectral properties of GUE and GOE and related ensembles 
on the local scale (see [0, |34| for the initial results for Gaussian ensembles and 
II f° r their generalizations). 

The case of band random matrices is different because the probability distri- 
bution of the ensemble H {n ^ (2.4) is no more invariant under transformations 
of the coordinates. One of the possible ways to study the spectral properties of 
jj(n,b) j g £ f n ow th e resolvent expansions approach well-known in theoretical 
physics (see, for exampl e A rigorous version of it has been developed in a 



series of papers |21| |22| g0[ . 

In frameworks of the resolvent approach (see |2(]] for details), one con- 
siders the correlation function C ny b (zi, z 2 ), Im Zj 7^ (3.4) in the limit when 
the dimension of the matrix N infinitely increases. Asymptotic expression for 
5(21,22) regarded in the limit z% = Ai + iO, z 2 = A2 — iO supplies one with 
the information about the local properties of eigenvalue distribution provided 
Ai — A2 = 0(N^ 1 ). Indeed, according to (2.7), the formal definition of the 
eigenvalue density p n .b(X) = <j' n b (X) is 

Pn,b{\) = j { [/„,&(A + iO) - /„, b (A - iO)] . 
Then one can consider expression 

fln,6(Ai,A 2 ) = -- S 1 S 2 C nib {X 1 +iS 1 0,\ 2 + iS 2 0) 

S 1 ,S 2 --l,+l 

as the correlation function of p n ,b- In general, even if R nt b can be rigorously 
determined, it is difficult to carry out the direct study of it. Taking into account 
relation (3.5), one can pass to more simple expression 

S n , 6 (A 1 ,A 2 ) = -j^ J2 5 1 S 2 S{X 1 +i5 1 0,\ 2 + iS 2 0) (6.1) 

5 U S 2 = -1, + 1 

and assume that it corresponds to the leading term of i? n ,f>(Ai, A2) in the limit 
n, b — > 00. 

In present section we follow the same heuristic scheme. It should be noted 
that for Wigner random matrices this approach is justified by the study of the 
simultaneous limiting transition N — -> oo,ImZj — > in the studies of C/v (21,22) 
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Theorem 6.1. 

Let S(zi,z 2 ) is given by (3.6). Assume that function uf(p) is such that there 
exist positive constants c\ , <5 and v > 1 that 

u F (p)=u F (0)-c 1 \p\ v +o(\p\") (6.2) 

for all p such that \p\ < 5, S — > 0. Then 

^,^ = ^^-^,(1 + 0(1)) (6.3) 

for Xj,j = 1,2 satisfying 

Ai,A 2 -» A e (-2y/v,2y/v). (6.4) 



Proof o/ Theorem 6.1. 

Let us start with the terms of (6.1) that correspond to 6\S 2 = — 1- It follows 
from (2.9) that 

\ — vw\W 2 = . (6.5) 

W\ — V02 

Also for the real and imaginary parts of w(X + iO) = r(A) + ip(A), we have 

2 A 2 2 4w-A 2 

(here and below we omit variable A). This implies existence of the limits w{z\) = 
w{z 2 ) for (6.4). One can easily deduce from (6.5) that in the limit (6.4) 

1 - vw( Zl )w(z 2 ) = X ^ lp ^ = ( 6 - 7 ) 

Also we have that 

(1 - vw 2 )(l - vw 2 2 ) = 2 - 2w(t 2 - p 2 ) = Avp 2 . (6.8) 
Now let us consider Q(z\, z 2 ) (3.8) and write that 

W 2 W 2 Up(p) 



Q(zi,z 2 ) = ±-( f + f ) 



2<ip = Ii + I 2 . 
[1 - vwiw 2 uf{p)] 



Relations (6.5) and (6.7) imply equality (cf. (3.9)) 

[1 - vw lW2 u F {p)] 2 = [u F (p) - 1] 2 (1 + o(l)). (6.9) 
Since u(t) is monotone, then 

liminf p£R \ ( _ 5i5) [u F (p) - l] 2 > 0. 
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This means that I 2 < 00 in the limit (6.4). 

Regarding (6.7), we can write that in the limit (6.4) 



Ii 



(27r) 1 w\w\uf{p) 



s (1 — vw\w 2 + VW1W2 [uf(p) — 1]) 
Then we derive relation 



7 dp 



(2Trv)- 1 u F { P )(l + o(l)) 

(£E£ + Mp) - 1] 



dp. 



Ii(Ai + iO, A 2 - iO) + Ii(Ai - iO, A 2 + iO) = 



■5 [Up 
8 



u F {p) - If 



2p 



r u F (p)(l + o(l))dp, 



(6.10) 



where o(l) corresponds to (6.9) regarded in the limit (6.4). 
Now let us use condition (6.2) and observe that 



5 c 2 p 2 » 



o{p 2 



D 2 



8 [c 2 p 2l/ + o(p 2u ) + D 2 } 



;dp ■ 



ttD 2 - 1 /" J 



8D- 



c\s 2v + o{s 2v ) - 1 
[c 2 s 2 " + o{s 2v ) + lf 



ds, 



where we denoted D = \\\ — A 2 |/(2p) and o{p 2u ) corresponds to the limit 5^0 
(6.2). Now it is clear that if we take 5 such that <5|Ai — A 2 | _1,/!/ — ► 00, we obtain 
asymptotically 

(2o) 2 ~ 1 ^ 
h+h=AB v { Cl )- - r 



|Ai-A 2 | 2 -V- 



where 



B v [ci) = 



1 



2ttc 



ds 



1 



- 2 



ds 



(1 + s 2 -) 2 



(6.11) 
(6.12) 



To prove relation (6.3), it remains to consider the sum 

J(Ai + iO, A 2 + iO) + J(Ai - iO, A 2 - iO). 

It is easy to observe that relations of the form (6.8) imply boundcdness of this 
sum in the limit (6.4) 

Now gathering relations (6.8) and (6.11), we derive that 



S,,(A„A 2 ), = i ^p— + „(!)) 



(6.13) 



This proves (6.3).D 



Let us discuss two consequences of Theorem 6.1. Let us assume first that 
u(t) is such that 

J t 2 u{i)At < 00. (6.14) 



u 2 
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Then (6.2) holds with v = 2 and ci = Regarding the right-hand side of (3.5) 
in the limit (6.4) with Xj = A + rj/N, j = 1, 2, we obtain asymptotic relation 



where 



B 2 (u 2 ) = f = -^=r (*) r (T) . (6.156) 



Now let us assume that (6.14) is not true. Suppose that there exists such 
1 < v' < 2 that 

u{t) = 0{\t\^ v ') as t^oo. (6.16) 

Then one can easily derive that (6.2) holds with v = v' . This follows from 
elementary computations based on equalities 



up(p) = u_f(0) — / (1 — cospt) u{t)&t 

J — oo 

and 

- r (l-co S y)u(yp- 1 )dy=\pf f°° {1 ~™ V) dy + o(\pf), p - 0. 
P J-oo J-oo |y| 

Therefore, if (6.16) holds, then 

The form of asymptotic expressions (6.15a) and (6.17) coincides with that 
determined by Altshuler and Shklovski for the spectral correlation function of 
band random matrices (see p7f for this and similar results). In these works, the 
factor | ri — r2|~ 3 / 2 appeared instead of usual for random matrices expression 
\t\ — r2\~ 2 (see (3.10)). This has been interpreted as the evidence of (relatively) 
localized eigenvectors of H^ n ' b ' in the limit 1 <C b <C n with the localization 
length b 2 /n. Let us note that the asymptotic expressions similar to (6.15) have 
also appeared in the recent work |33| , where the band random matrix ensemble 
ff( n , b ) was considered under condition (6.14). However it should be stressed 
that no explicit expressions like (3.6) and (6.15) were obtained neither in p7fl 
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7 Summary 

We consider a family of random matrix ensembles of the band-type 

form. More precisely, we are related with real symmetric N x N matrices, 
N = 2n + 1, whose entries are jointly independent Gaussian random variables 
with zero mean value. The band-type form means that the variance of the 
matrix entries H^ n ' b '(x, y) is proportional to u (^y^) > 0. 
We study asymptotic behavior of the correlation function 

C n ,b(zi,Z 2 ) = ~E>fn,b{ Z l)fn,b(z 2 ) - ~Ef n>b (zi )E/ n>6 (zi ) , 

where f n ,b{z) is the normalized trace (G^ n ' b '(z)') of the resolvent of H^ n ' b \ 
We have proved that if Im zj is large enough, then in the limit 1 <C b <C N 1 / 3 

Cn,b{zi,Z 2 ) = J^S{Z1, Z 2 ) + O 

We have found explicit form of the leading term S{z\ 1 z 2 ) in this limit. Assuming 
that expression £ n) (,(Ai, \ 2 ) (6.1) is closely related with the correlation function 
of the eigenvalue density, we have studied it in the limit N, b — > oo and Ai — X 2 = 
(n-r 2 )/N. 

Our main conclusion is that the limiting expression for S n! b exhibits different 
behavior depending on the rate of decay of u(t) at infinity. 
If J t 2 u(t)dt < oo, then (6.1) is given by 



-^F^F (1 + O(1)) ' c>0 - 

If u{t) = 0(|i j _I/ ) with 1 < v < 2, then the asymptotic expression for (7.1) is 
proportional to 

n 1 - 1 /- 1 



b Iri-ral 2 - 1 ^' 

In both cases the exponents do not depend on the particular form of the function 
u(t). Moreover, in the first case the exponents do not depend on u at all. This 
can be regarded as a kind of spectral universality for band random matrix 
ensembles. On can conject that these characteristics also do not depend on the 
probability distribution of the random variables a(x, y) (2.1). 

Our results show that S(z\, z 2 ) determines at least two scales of universality 
in the local spectral properties of band-type random matrices. These scales 
coincide with those detected in theoretical physics for the (relative) localization 
length and density-density correlation function for these ensembles f27j| . In the 
papers cited also the third scale when u(i) = 0(|£| -7 ) with 7 £ (1/2, 1) has been 
observed. It been shown to produce the asymptotics 7V~ 2 |ri — r 2 \~ 2 which is 
typical for "full" random matrices like GOE ||, [l4|, Unfortunately, this 
asymptotic regime for band random matrices is out of reach of our technique. 
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